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Results on two-particle angular correlations in jet cones and on multiplicity fluc- 
*Vj ' tuations in one- and two- dimensional angular intervals, delivered by three exper- 

iments (DELPHI, L3 and ZEUS, at y/s from few to 183 GeV) are compared to 
present existing analytical QCD calculations, using the LPHD hypothesis. Two 
different types of functions have been tested. While the differentially normalized 
correlation functions show substantial deviations from the predictions, a globally 
normalized correlation function agrees surprisingly well. The role of the QCD pa- 
rameters as, A and Uf is discussed. The necessity to include full energy- momentum 
conservation into the analytical calculations is stressed. 
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1 Introduction 



Phenonienological models which try to describe multihadron production in 

fpH, e+e~ reactions have to cut off the parton cascade at some scale Qo > 1 GeV. 

O ,' The following non-perturbative hadronisation phase is "handled" with many 

D , parameters such that the connection to QCD is "diluted". It has been sug- 

"^ ' gested to extend instead the parton evolution calculated using perturbative 

QCD down to a lower mass scale. The multihadron final states are then com- 
pared directly with the multiparton final statesLU by using the concept of Local 
Parton Hadron Duality (LPHD) □. The obvious difficulty is the overlap of 
5^ I perturbative and non-perturbative QCD. It is unclear what one could expect 

when entering a region in which nobody knows what to apply. Therefore an 
experimental input is vital and the only way to learn what on really can ap- 
ply. The main theoretical effort for evaluating multiparton correlations in the 
framework of QCD has been based on the Double Log Approximation (DLA) 
liJo. Comparisons with experimental data have to cope with substantial sim- 
plifications in the calculations of the perturbative part which are justified only 
at asymptotic energies, as well as with the question of how far the LPHD 
hypothesis is valid. 

Detailed prescriptions for multiparton angulaL correlations in cones us- 
ing DLA have been proposed by three groups QoH, which calculated parton 
correlations produced in gluon cascades radiated off the initial parton. The 



analytical predictions considered thereby involve only one adjustable param- 
eter, namely the QCD scale A. The probability for gluon bremsstrahlung in 
DLA reads as follows: 

M(fc)d3fc = c.7o4#^^ 7o^ = 6asA (1) 

where Cg — 1 and Cq — „, p and k are the 3-momcnta of the parent parton and 
the radiated gluon, 8^^ the angle of emission of the gluon, ^pk the azimuthal 
angle of the gluon around p and 70 is dependent on the pt of the gluon. The in- 
clusive n-particle densities pn{ki, ^2, fc„) [ki is the S-momentum of the i-th 

particle) are obtained by applying the generating functional technique lil which 
has been developed for QCD jets first byl3'EI. The calculations have been car- 
ried out in DLA where the integrals involved are performed only in phase space 
regions with dominant contributions given by the singularities of (1). Among 
the simplifications that had to be used are neglecting energy-momentum con- 
servation and qq production, and assuming well developed cascades at very 
high energies. Angular ordering, however, has been taken into account. 

The predictions ofQoIrl and the present study use the lowest order QCD 
relations between the coupling as and the QCD scale A. 

«^ = ^^_ /32 ^ i2(iin, - -n,)-i (2) 

6 ln(2) ^33^^ ^ ' 

The aim of this study is to compare experimental measurements of angu- 
lar correlations in cones around the jet axis as well as multiplicity fluctuations 
in one- and two- dimensional angular intervals with the available theoretical 
QCD predictions thus helping the theorists to assess how well their analyti- 
cal calculationisand what improvements could be done. First experimental 
measurementsH'cLill revealed substantial deviations, but also encouraging agree- 
ments, rnm rn 

This study uses resuhs of theLEP experiments DELPHlliH-BJ and L3EJ 
and the HERA experiment ZEUSLLj. Section 2 contains information about the 
experimental data used for the comparison. In sections 3 and 4 the theoretical 
framework is sketched and the analytical calculations are compared with the 
experimental measurement. Section 5 contains the final discussion and the 
summary. 

2 The data samples 

All three experimentaxonsklfired use^ecific cuts for hadronic events and track 
quality (DELPHlli2Hi3, L3li^, ZEUSLLs) Various corrections were applied using 



events generated from Monte Carlo simulations! Hnr'^H^. In the following 
only corrected data will be shown. 

DELPHI : The analysis uses about 600k e+e^ interactions (after cuts) at 
y/s = 91 GeV . A sample of about 1200 high energy events at y/s — 183 GeV 
collected in 1997 is used to investigate the energy dependence. 

L3 : About 1 million e+e" interactions at a/s = 91 GeV are used. 

The event axis is the sphericity axis both in DELPHI and L3. 

ZEUS : The data of this experiment {ep interactions at HERA) are 4 
samples with different ranges of Q^: 85k events (10 < Q'^ < 20 GeV), 28k 
events (Q^ > 100 GeV), 0.9k events (Q^ > loOO GeV) and 0.3k events with 
(Q^ > 2000 GeV). The last sample corresponds to y/s ~ 62 GeV, which comes 
near to LEPl energies. The event axis is defined using the initial quark mo- 
mentum of the 7*q collision in the quark-parton model. 

3 Multiplicity fluctuations in 1- and 2- dimensional ring regions 
around jet cones 

The theoretical calculations treat correlations between partons emitted within 
an angular window defined by two angles d and Q. The parton and particle 
density correlations (fluctuations) in this window are described by normalized 
factorial moments of order n: 

F<-'(e..). JnL^»>.;"'(n..^.,o.) p, 

J lL=iaS2feP^^H^^fc) 

The integrals extend over the window chosen. The angular windows considered 
here are either rings around the jet axis with mean opening angle and half 
width i} in the case of 1 dimension (D = 1), or cones with half opening angle d 
around a direction (0, $) with respect to the jet axis in the case of 2 dimensions 
{D = 2). At sufficiently large jet energies, the parton flow in these angular 
windows is dominated by parton avalaches caused by gluon bremsstrahlung off 
the initial quark. 

Ref. H derived their predictions explicitly for cumulant moments C^"' lij, 
whereas B and El obtained similar expressions for the factorial moments F^"^' . 

Ear the normalized cumulant moments C'"'cl and the factorial moments 
F^^^ElO, the following prediction has been made: 

C(")(e,i?) or F(")(e,?9) - f^'] " (4) 

All 3 referencesuO'EI give in the high energy limit and for large i} < Q the same 



linear approximation for the exponents </)„: 

<^„ « [n - 1)D - L - i j 70 (5) 

For fixed as (along the parton shower) eq. (5) is asymptotically valid for all 
angles ■&. 

_When the running of as with t? in the parton cascade is taken into account, 
ref. □ obtained 

0„ = (n-l)i?-27o(n-c^(e,n))/e (6) 

where D = 1 for ring regions and Z? = 2 for cones, 



1 
2n 
and 



uj{€,n)=7i^/Y~'e{l-^^\n{\~e)) (7) 



in(e/.) ^^^ 



ln(Pe/A) 

where P « -^ is the momentum of the initial parton. 

The dependence on the QCD parameters a^ or A enters in the above 
equations via 70 and e that are determined by the scale Q ~ PQ. In the 
present study it is, for Q = 25°, about 20 GeV for -/s=91.1 GeV and about 
38 GeV for ^ w 183 GeV. 

The corresponding predictions of refs.El (eq. 9) andEl (eq. 10) are analyti- 
cally different, but numerically similar: 

2 



270 n^-1 



0„ - {n -l)D~^ (1 - VT~e) (9) 

en 

0„ ^{n- l)D - ^70 (1 + 4^e) (10) 

These relations depend also on the number of flavours («/). Since equ. 2 
emerges only from "one loop" calculations, the parameter A is not the universal 
^MS ' ^^1 only an effective parameter Ae// . But also in this approximation a^ 
is running having a scale dependence l/ln((3^/A^). The running of a^ during 
the process of jet cascading is implicitly taken into account in (6), (9) and (10) 
by the dependence of 0„ on e (or §). In theory this causes a bending of F*^"^ 
when approaching smaller values of § (larger e) . 
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Figure 1: Factorial moments in l-dim rings are compared with three Monte Carlo models, 
both on hadronic and partonic level. Please note that the variable e (equ. 8) is here denoted 
by z. 



In fig. 1 (L3) the factorial moments of orders 2,3,4 and 5 are compared to 
various Monte Carlo simulations at hadronic and partonic level. For all orders 
n the hadronic levels agree very well with the data (this is generally the case), 
while the partonic level disagrees with that of the hadronic level, so that LPHD 
seems to be not valid in this case. This is somewhat contradicted by ZEUS 
lU, which claim an approximate validity for LPHD for Q^ > 100 GeV^. When 
normalizing these moments by the respective Jialue at e = the agreement 
between parton and hadron levels is improvedt!-^. It is pointed to the fact 
that the shape of the correlation fimctions on the partonic level depends very 
much on the cut-off parameter QqI— ' (in fig. 1 Qo=l GeV, a rather high value 
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Figure 2: Factorial moments_Ln_a) 1-dim rings and b) 2-dim side cones are compared with 
analytical calculations of refsQUlHI, eqs. (6) (solid lines), (9) (dashed lines), and (10) (dotted 
lines), for a cone opening angle of © = 25°. As a consistency test, 1-dim and 2-dim factorial 
moments are compared with same values of QCD parameters: note the different vertical 
scales. The values for orders 2 to 5 are indicated in all figures. The statistical errors are 
shown by the error bars, the shaded regions indicate the systematic errors. 



compared to the original understanding of LPHDa). 

Fig. 2 (DELPHI) shows the normahzcd factorial moments of orders 2, 3, 4 
and 5 together with the predictions of refs.Qi3'cl, in one- and two- dimensional 
angular intervals (i.e. rings and side cones) for A— 0.15 GeV and n/=3. They 
are not described well by the theoretical predictionsa\B'a 

For the 1-dim case (fig. la) the predictions lie below the data for not 
too large e, differing also in shape. Choosing nj = 5 the discrepances will 
increase liiJ, choosing smaller values of A (e.g. = 0.04 GeV) will reduce the 
discrepancies for small e Iiil^li2l, especially for lower orders n. The angular 
correlations in 2 dimensions (fig. 2b) show a different (disagreeing) behaviour 
for the lower order moments n < 4, where the predictions lie above the data. 
The higher moments F^^'' and F'-^'> have similar features in the 1-dim and 2- 
dim case. In both cases the data lie above the predictions at small e and bend 
below them at larger e. 

Similar conclusions for y^=91 GeV, for the 1-dim ring regions, are drawn 
by the L3 collaboration. 

It is not possible to find one set of QCD parameters A and Uf which simul- 
taneously minimize the discrepancies between data and predictions for moments 
of all orders 2,3,4 (^""^d 5 in both the 1-dim. and 2-dim. cases. 

Fig. 3 (DELPHI) shows a comparison with high energy data at \/s=183 GeV 
(with a mean energy of y^=175 GeV and the corresponding predictions ac- 
cording to equ. (6), where the energy dependence enters via the parameter 
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Figure 3: The energy dependence of 

the normaUzed factorial moment of order 
2 y/s=91.1 GeV: data (open circles), in 
full agreement with JETSET, and predic- 
tionQ dashed line and Y^=pL75 GeV: data 
(full circles) and predictionQ solid line, for 
n/ = 3 and A=0.04 GeV. The full trian- 
gles denote the high energy JETSET sim- 
ulation. 



Figure 4: Values of 7o obtained from 
fitting the linear approximation eq. (5) vs. 
§ the 1-dim case (open circles) and the 
2-dim case (full circles), for orders n = 
2, 3, 4, 5. The measured values of 7q are 
also compared with those obtained from 
JETSET at generator level: open triangles 
(1-dim) and full triangles (2-dim). 



70- For larger values of e there is better agreement, the statistical errors of 
the high energy data, however, are substantial. The relative increase of the 
predicted moments agrees qualitatively with that of the JETSET model. 

The QCD parameter 70 is discussed in fig. 4 (DELPHI), where the numer- 
ical values of 7o derived from the measured slopes (f>n (by fitting equ. (5) 
for e < 0.1) are given for the orders n = 2, 3, 4, 5. From the present theoretical 
understanding, 70 is expected to be independent of n. This is indicated by 
horizontal lines in fig. 4. It has to be pointed out that the average measured 
values of 7o s-re not too far from the expectation. The n- dependence ob- 
served, however, is not described by the calculations. The measured values of 
7q agree, however, extremely well with the corresponding values obtained 
from JETSET, as can be seen in fig. 4. 



4 Predictions on 2-particle angular correlations in jet cones 

Theoretical predictionsu concerning the emission of two partons with a relative 
angle ■(?i2 - within a cone with half opening angle Q around the jet axis - have 
been evaluated using two correlation functions defined as follows : 



r(^i 



P2(^l 



Pi «)Pl(^12) 



r(^i2) 



P2{^ 



12 J 



712(6) 



(11) 



(12) 



with the correlation integralsEyp2 (^12) = /q (i'^fcid'^fc2P2(fci, fe)<5(i?i2— ^(fci, fe)) 
and 

Pi'^pi{-&12) = JQ(Pki(Pk2Pi{ki)pi{k2)5{di2-'d{ki,k2)) where pi(A:) is the sin- 
gle particle distribution and n(8) is the mean multiplicity of partons emitted 
into the 0-cone. The quantities in eqns. (11) and (12) exhibit very differ- 
ent structures. ^2(^12) consists of 2 terms ^2(^12) = C'2(i9i2) + Pi® Pi (^12) 
where only C2('!?i2) describes the genuine correlations, pi® Pl{'&l2)^ on the 
other side, is obtained from the single particle spectra. Consequently r(^i2) 
is given essentially by the normalized C2-term whereas it turns out that the 
dominant term of r(i?i2) is given by pi ® pi(i?i2). 

Distinct predictions for r(??i2) and f(^i2) have been evaluated which de- 
pend essentially only on the QCD parameters A and n/ Q. 

At high energy and for sufficiently large angles ^12 < 6 the following 
power law is expected: 

Q \ 0-570 



ri^.2) = ^^j (13) 

and the scale determining 70 is again given by Q « PQ (see section 3). 
For asymptotically high energies the quantity 

ln(r(7?i2)) 



In^ 



with 

In^ 



In-^ 



2/3(w(e, 2) - 2VT^) (14) 



/?2 = 12(11 --n/)-i (15) 



and w given by equ. 7, is expected to be independent of the cone opening angle 
Q and primary momentum P, meaning that it is a scaling function. 

The expected scaling properties of the quantity "^''^ ^'^'' (eqn. (14)) are 

tested in figs. 5 and 6. 

The dependence on the cone opening angle 8 is shown in figs. 5b (ZEUS) 
and 6a (DELPHI) (left side of figs 5, 6). It can be seen that for 15° < 6 < 60° 
the dependence on 8 is very weak already at ^/s = 91 GeV , in agreement with 
the predictions of eqn. (14). For small energies (fig 5b) there is no scaling for 
e < 0.3. The shape predicted by eqn. (14) differs appreciably from the measure- 
ment at LEP (fig. 6). The shape of the data is only similar to that predicted 
in the sense that it is rising and levelling off; the data are much smaller and 
fiatter. There is a "hook" in data at small e which is a refiection of momen- 
tum conservation. It is almost absent, however, in the ZEUS data. Thus at 
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Figure 5: The r2 correlation function (14) dependence on li) cone size © and a) 
lower energy range - and comparison with QCD predictionsQ. 
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Figure 6: a) The.rescaled function r2 (14) is plotted against the scaling variable e in order 
to test predictionQ, dashed line, b) The corrected data for = 45°, at 91 GeV (open circles) 
and at 183 GeV (full circles), using A = 0.15 GeV, are shown together with JETSET Monte 
Carlo calculations (open resp. full triangles). 
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f(e)=^i2f(^9i2)ln — 


ln(f(e)/6) 
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LEP energies the analytic QCD calculations do not describe quantitatively 
the 2-particle angular correlations r('i9i2). These differences are much less pro- 
nounced for the low energies (solid lines in fig. 5), but this could be also due 
to the different way to define the jet axis. The dependence of (14) on energy is 
shown in figs. 5a (ZEUS) and 6b (DELPHI) for a fixed value of 8 (right hand 
sides of figs. 5^6). The ZEUS data indicate an energy dependence only at low- 
est energies 1^3. This is in disagreement with the the observations of DELPHI 
which observes a strong energy dependence also at the highest energies. 

In equs. (16), (17) a scaling function y(e) is defined, which is predicted 
to be independent of and the primary momentum P. 

(16) 

ne) = -^^^7^=^^=^ = 2/3(l-0.5u;(e,2)) , 6 = 2/3Wln(l-^) (17) 

The expected scaling properties of Y{e) are tested in figs. 7 (ZEUS) and 
8 (DELPHI). From ^ ss 18 GeV upwards the distributions agree very well 
with the prediction in the region e > 0.2, therefore exhibiting scaling both 
in angle Q and energy. There is also good agreement with the coresponding 
JETSET simulations on both the partonic and hadronic level, which supports 
parton hadron duality - for the special function Y{e). Note that no arbitrary 
normalization has been applied in the above figures. 

Similar to the case of factorial moments (fig. 4), values of Jq-' have been 
extracted from r(i?i2) by fitting the DELPHI data to equ. (13) El. The 
corresponding values of as depending on the cone opening angle O are shown 
in fig. 9a. Similar to the case of factorial moments the data agree better with 
lower values of A. This is at variance with the function Y{e) where the data 
agree best with the value A w 0.3 GeV (fig. 9b). It has to be noted that DLA 
takes only the leading singularities in all cases considered which cauld lead to 
different redefinitions of the effective QCD parameters (e.g. Ae//)li^. 

5 Summary and conclusions 

Present available QCD predictions on angular correlations, based on the DLA 
approximation, with first order relationship between as and A, are checked, 
comparing them to relevant results of three collaborations (DELPHI, L3 and 
ZEUS). The experimental measurements, which agree generally well with Monte 
Carlo simulations on hadronic level and not so well with the partonic level cover 
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Figure 7: The function Y(e) (equ. 17) evolving with b) Q^ and a) cone opening angle 0. 
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Figure 9: 



2-particle angular correlations in cones and l-dim and 2-dim multiplicity fluc- 
tuations in angular intervals. It turns out that functions which contain mainly 
single particle terms {Y{e)) are predicted well in every respect. Also the extrac- 
tion of QCD parameters (i.e. ^eff) out of these measurements leads to reason- 
able results. Differentially normalised correlations (r(i5i2), F(§) that contain 
2 or more particles, however, show often strong disagreements in shapes. The 
basic qualitative features, on the other hand, are fulfilled. Here smaller values 
of Ae// (resp. 7o or a'^-^-^) are favoured as well as nf=3 rather than nf=5. 
Up to now no set of QCD parameters will minimise simultaneously the various 
discrepancies. 

The main shortcoming of the analytical calculations is the missing energy- 
momentmnixinservation adopted. Some "partial" introduction of NLO cor- 
rections ILyiij did not improve matters significantly. 

Considering, however, that the analytical calculation does not use any free 
parameter besides A or a^, that they treat all kinds of correlations functions 
in an universal manner, that the asymptotic energies might be still far away 
for some subtle functions and that, nevertheless, the basic features of the pre- 
dictions are always seen in the experimental measurement, further theoretical 



efforts are encouraged. 

We thank S.Chckanov, J.Fuster, W.Kittel and W.Ochs for valuable advice. 
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